Comprehensive Characterization of LAT1 Cholesterol-Binding Sites.

The human L-type amino acid transporter 1 (LAT1; SLC7A5), is an amino acid exchanger protein, primarily found in the blood-brain barrier, placenta, and testis, where it plays a key role in amino acid homeostasis. Cholesterol is an essential lipid that has been highlighted to play a role in regulating the activity of membrane transporters, such as LAT1, yet little is known about the molecular mechanisms driving this phenomenon. Here we perform a comprehensive computational analysis to investigate cholesterol's role in LAT1 structure and function, focusing on four cholesterol-binding sites (CHOL1-4) identified in a recent LAT1-apo inward-open conformation cryo-EM structure. Through a series of independent molecular dynamics (MD) simulations, molecular docking, MM/GBSA free energy calculations, and other analysis tools, we explored the interactions between LAT1 and cholesterol. Our findings suggest that CHOL3 forms the most stable and favorable interactions with LAT1. Principal component analysis (PCA) and center of mass (COM) distance assessments show that CHOL3 binding stabilizes the inward-open state of LAT1 by preserving the spatial arrangement of the hash and bundle domains. Additionally, we propose an alternative cholesterol-binding site for originally assigned CHOL1. Overall, this study improves the understanding of cholesterol's modulatory effect on LAT1 and proposes candidate sites for the discovery of future allosteric ligands with rational design.

Development of a human genetics-guided priority score for 19,365 genes and 399 drug indications.

Studies have shown that drug targets with human genetic support are more likely to succeed in clinical trials. Hence, a tool integrating genetic evidence to prioritize drug target genes is beneficial for drug discovery. We built a genetic priority score (GPS) by integrating eight genetic features with drug indications from the Open Targets and SIDER databases. The top 0.83%, 0.28% and 0.19% of the GPS conferred a 5.3-, 9.9- and 11.0-fold increased effect of having an indication, respectively. In addition, we observed that targets in the top 0.28% of the score were 1.7-, 3.7- and 8.8-fold more likely to advance from phase I to phases II, III and IV, respectively. Complementary to the GPS, we incorporated the direction of genetic effect and drug mechanism into a directional version of the score called the GPS with direction of effect. We applied our method to 19,365 protein-coding genes and 399 drug indications and made all results available through a web portal.

Genome-wide prediction of pathogenic gain- and loss-of-function variants from ensemble learning of a diverse feature set.

Gain-of-function (GOF) variants give rise to increased/novel protein functions whereas loss-of-function (LOF) variants lead to diminished protein function. Experimental approaches for identifying GOF and LOF are generally slow and costly, whilst available computational methods have not been optimized to discriminate between GOF and LOF variants. We have developed LoGoFunc, a machine learning method for predicting pathogenic GOF, pathogenic LOF, and neutral genetic variants, trained on a broad range of gene-, protein-, and variant-level features describing diverse biological characteristics. LoGoFunc outperforms other tools trained solely to predict pathogenicity for identifying pathogenic GOF and LOF variants and is available at https://itanlab.shinyapps.io/goflof/ .

Substrate binding and inhibition of the anion exchanger 1 transporter.

Anion exchanger 1 (AE1), a member of the solute carrier (SLC) family, is the primary bicarbonate transporter in erythrocytes, regulating pH levels and CO2 transport between lungs and tissues. Previous studies characterized its role in erythrocyte structure and provided insight into transport regulation. However, key questions remain regarding substrate binding and transport, mechanisms of drug inhibition and modulation by membrane components. Here we present seven cryo-EM structures in apo, bicarbonate-bound and inhibitor-bound states. These, combined with uptake and computational studies, reveal important molecular features of substrate recognition and transport, and illuminate sterol binding sites, to elucidate distinct inhibitory mechanisms of research chemicals and prescription drugs. We further probe the substrate binding site via structure-based ligand screening, identifying an AE1 inhibitor. Together, our findings provide insight into mechanisms of solute carrier transport and inhibition.

Exploring the Druggable Conformational Space of Protein Kinases Using AI-Generated Structures.

Protein kinase function and interactions with drugs are controlled in part by the movement of the DFG and ɑC-Helix motifs, which enable kinases to adopt various conformational states. Small molecule ligands elicit therapeutic effects with distinct selectivity profiles and residence times that often depend on the kinase conformation(s) they bind. However, the limited availability of experimentally determined structural data for kinases in inactive states restricts drug discovery efforts for this major protein family. Modern AI-based structural modeling methods hold potential for exploring the previously experimentally uncharted druggable conformational space for kinases. Here, we first evaluated the currently explored conformational space of kinases in the PDB and models generated by AlphaFold2 (AF2) (1) and ESMFold (2), two prominent AI-based structure prediction methods. We then investigated AF2's ability to predict kinase structures in different conformations at various multiple sequence alignment (MSA) depths, based on this parameter's ability to explore conformational diversity. Our results showed a bias within the PDB and predicted structural models generated by AF2 and ESMFold toward structures of kinases in the active state over alternative conformations, particularly those conformations controlled by the DFG motif. Finally, we demonstrate that predicting kinase structures using AF2 at lower MSA depths allows the exploration of the space of these alternative conformations, including identifying previously unobserved conformations for 398 kinases. The results of our analysis of structural modeling by AF2 create a new avenue for the pursuit of new therapeutic agents against a notoriously difficult-to-target family of proteins. Significance Statement: Greater abundance of kinase structural data in inactive conformations, currently lacking in structural databases, would improve our understanding of how protein kinases function and expand drug discovery and development for this family of therapeutic targets. Modern approaches utilizing artificial intelligence and machine learning have potential for efficiently capturing novel protein conformations. We provide evidence for a bias within AlphaFold2 and ESMFold to predict structures of kinases in their active states, similar to their overrepresentation in the PDB. We show that lowering the AlphaFold2 algorithm's multiple sequence alignment depth can help explore kinase conformational space more broadly. It can also enable the prediction of hundreds of kinase structures in novel conformations, many of whose models are likely viable for drug discovery.

Targeting SLC transporters: small molecules as modulators and therapeutic opportunities.

Solute carrier (SLCs) transporters mediate the transport of a broad range of solutes across biological membranes. Dysregulation of SLCs has been associated with various pathologies, including metabolic and neurological disorders, as well as cancer and rare diseases. SLCs are therefore emerging as key targets for therapeutic intervention with several recently approved drugs targeting these proteins. Unlocking this large and complex group of proteins is essential to identifying unknown SLC targets and developing next-generation SLC therapeutics. Recent progress in experimental and computational techniques has significantly advanced SLC research, including drug discovery. Here, we review emerging topics in therapeutic discovery of SLCs, focusing on state-of-the-art approaches in structural, chemical, and computational biology, and discuss current challenges in transporter drug discovery.

Conserved allosteric inhibition mechanism in SLC1 transporters.

Excitatory amino acid transporter 1 (EAAT1) is a glutamate transporter belonging to the SLC1 family of solute carriers. It plays a key role in the regulation of the extracellular glutamate concentration in the mammalian brain. The structure of EAAT1 was determined in complex with UCPH-101, apotent, non-competitive inhibitor of EAAT1. Alanine serine cysteine transporter 2 (ASCT2) is a neutral amino acid transporter, which regulates pools of amino acids such as glutamine between intracellular and extracellular compartments . ASCT2 also belongs to the SLC1 family and shares 58% sequence similarity with EAAT1. However, allosteric modulation of ASCT2 via non-competitive inhibitors is unknown. Here, we explore the UCPH-101 inhibitory mechanisms of EAAT1 and ASCT2 by using rapid kinetic experiments. Our results show that UCPH-101 slows substrate translocation rather than substrate or Na+ binding, confirming a non-competitive inhibitory mechanism, but only partially inhibits wild-type ASCT2. Guided by computational modeling using ligand docking and molecular dynamics simulations, we selected two residues involved in UCPH-101/EAAT1 interaction, which were mutated in ASCT2 (F136Y, I237M, F136Y/I237M) in the corresponding positions. We show that in the F136Y/I237M double-mutant transporter, 100% of the inhibitory effect of UCPH-101 could be restored, and the apparent affinity was increased (Ki = 4.3 µM), much closer to the EAAT1 value of 0.6 µM. Finally, we identify a novel non-competitive ASCT2 inhibitor, through virtual screening and experimental testing against the allosteric site, further supporting its localization. Together, these data indicate that the mechanism of allosteric modulation is conserved between EAAT1 and ASCT2. Due to the difference in binding site residues between ASCT2 and EAAT1, these results raise the possibility that more potent, and potentially selective ASCT2 allosteric inhibitors can be designed .

Structural Analysis of Genomic and Proteomic Signatures Reveal Dynamic Expression of Intrinsically Disordered Regions in Breast Cancer and Tissue.

Structural features of proteins capture underlying information about protein evolution and function, which enhances the analysis of proteomic and transcriptomic data. Here we develop Structural Analysis of Gene and protein Expression Signatures (SAGES), a method that describes expression data using features calculated from sequence-based prediction methods and 3D structural models. We used SAGES, along with machine learning, to characterize tissues from healthy individuals and those with breast cancer. We analyzed gene expression data from 23 breast cancer patients and genetic mutation data from the COSMIC database as well as 17 breast tumor protein expression profiles. We identified prominent expression of intrinsically disordered regions in breast cancer proteins as well as relationships between drug perturbation signatures and breast cancer disease signatures. Our results suggest that SAGES is generally applicable to describe diverse biological phenomena including disease states and drug effects.

Alanine serine cysteine transporter (ASCT) substrate binding site properties probed with hydroxyhomoserine esters.

The glutamine transporter ASCT2 is highly overexpressed in cancer cells. Block of glutamine uptake by ASCT2 is a potential strategy to inhibit growth of cancer cells. However, pharmacology of the ASCT2 binding site is not well established. In this work, we report the computational docking to the binding site, and the synthesis of a new class of ASCT2 inhibitors based on the novel L-hydroxyhomoserine scaffold. While these compounds inhibit the ASCT2 leak anion conductance, as expected for competitive inhibitors, they did not block leak conductance in glutamate transporters (EAAT1-3 and EAAT5). They were also ineffective with respect to subtype ASCT1, which has >57% amino acid sequence similarity to ASCT2. Molecular docking studies agree very well with the experimental results and suggest specific polar interactions in the ASCT2 binding site. Our findings add to the repertoire of ASCT2 inhibitors and will aid in further studies of ASCT2 pharmacology.

Functional genomics of OCTN2 variants informs protein-specific variant effect predictor for Carnitine Transporter Deficiency.

Genetic variants in SLC22A5, encoding the membrane carnitine transporter OCTN2, cause the rare metabolic disorder Carnitine Transporter Deficiency (CTD). CTD is potentially lethal but actionable if detected early, with confirmatory diagnosis involving sequencing of SLC22A5. Interpretation of missense variants of uncertain significance (VUSs) is a major challenge. In this study, we sought to characterize the largest set to date (n = 150) of OCTN2 variants identified in diverse ancestral populations, with the goals of furthering our understanding of the mechanisms leading to OCTN2 loss-of-function (LOF) and creating a protein-specific variant effect prediction model for OCTN2 function. Uptake assays with 14C-carnitine revealed that 105 variants (70%) significantly reduced transport of carnitine compared to wild-type OCTN2, and 37 variants (25%) severely reduced function to less than 20%. All ancestral populations harbored LOF variants; 62% of green fluorescent protein (GFP)-tagged variants impaired OCTN2 localization to the plasma membrane of human embryonic kidney (HEK293T) cells, and subcellular localization significantly associated with function, revealing a major LOF mechanism of interest for CTD. With these data, we trained a model to classify variants as functional (>20% function) or LOF (<20% function). Our model outperformed existing state-of-the-art methods as evaluated by multiple performance metrics, with mean area under the receiver operating characteristic curve (AUROC) of 0.895 ± 0.025. In summary, in this study we generated a rich dataset of OCTN2 variant function and localization, revealed important disease-causing mechanisms, and improved upon machine learning-based prediction of OCTN2 variant function to aid in variant interpretation in the diagnosis and treatment of CTD.

Describing inhibitor specificity for the amino acid transporter LAT1 from metainference simulations.

The human L-type amino acid transporter 1 (LAT1; SLC7A5) is a membrane transporter of amino acids, thyroid hormones, and drugs such as the Parkinson's disease drug levodopa (L-Dopa). LAT1 is found in the blood-brain barrier, testis, bone marrow, and placenta, and its dysregulation has been associated with various neurological diseases, such as autism and epilepsy, as well as cancer. In this study, we combine metainference molecular dynamics simulations, molecular docking, and experimental testing, to characterize LAT1-inhibitor interactions. We first conducted a series of molecular docking experiments to identify the most relevant interactions between LAT1's substrate-binding site and ligands, including both inhibitors and substrates. We then performed metainference molecular dynamics simulations using cryoelectron microscopy structures in different conformations of LAT1 with the electron density map as a spatial restraint, to explore the inherent heterogeneity in the structures. We analyzed the LAT1 substrate-binding site to map important LAT1-ligand interactions as well as newly described druggable pockets. Finally, this analysis guided the discovery of previously unknown LAT1 ligands using virtual screening and cellular uptake experiments. Our results improve our understanding of LAT1-inhibitor recognition, providing a framework for rational design of future lead compounds targeting this key drug target.

Structural signatures: a web server for exploring a database of and generating protein structural features from human cell lines and tissues.

Structural features of proteins provide powerful insights into biological function and similarity. Specifically, previous work has demonstrated that structural features of tissue and drug-treated cell line samples can be used to predict tissue type and characterize drug relationships, respectively. We have developed structural signatures, a web server for annotating and analyzing protein features from gene sets that are often found in transcriptomic and proteomic data. This platform provides access to a structural feature database derived from normal and disease human tissue samples. We show how analysis using this database can shed light on the relationship between states of single-cell RNA-sequencing lung cancer samples. These various structural feature signatures can be visualized on the server itself or downloaded for additional analysis. The structural signatures server tool is freely available at https://structural-server.kinametrix.com/.

Exploring the conformational diversity of proteins.

An artificial intelligence-based method can predict distinct conformational states of membrane transporters and receptors.

Identification of discriminative gene-level and protein-level features associated with pathogenic gain-of-function and loss-of-function variants.

Identifying whether a given genetic mutation results in a gene product with increased (gain-of-function; GOF) or diminished (loss-of-function; LOF) activity is an important step toward understanding disease mechanisms because they may result in markedly different clinical phenotypes. Here, we generated an extensive database of documented germline GOF and LOF pathogenic variants by employing natural language processing (NLP) on the available abstracts in the Human Gene Mutation Database. We then investigated various gene- and protein-level features of GOF and LOF variants and applied machine learning and statistical analyses to identify discriminative features. We found that GOF variants were enriched in essential genes, for autosomal-dominant inheritance, and in protein binding and interaction domains, whereas LOF variants were enriched in singleton genes, for protein-truncating variants, and in protein core regions. We developed a user-friendly web-based interface that enables the extraction of selected subsets from the GOF/LOF database by a broad set of annotated features and downloading of up-to-date versions. These results improve our understanding of how variants affect gene/protein function and may ultimately guide future treatment options.

Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: The Multi-Targeting Drug DREAM Challenge.

A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ('polypharmacology'). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Rational design of ASCT2 inhibitors using an integrated experimental-computational approach.

ASCT2 (SLC1A5) is a sodium-dependent neutral amino acid transporter that controls amino acid homeostasis in peripheral tissues. In cancer, ASCT2 is up-regulated where it modulates intracellular glutamine levels, fueling cell proliferation. Nutrient deprivation via ASCT2 inhibition provides a potential strategy for cancer therapy. Here, we rationally designed stereospecific inhibitors exploiting specific subpockets in the substrate binding site using computational modeling and cryo-electron microscopy (cryo-EM). The final structures combined with molecular dynamics simulations reveal multiple pharmacologically relevant conformations in the ASCT2 binding site as well as a previously unknown mechanism of stereospecific inhibition. Furthermore, this integrated analysis guided the design of a series of unique ASCT2 inhibitors. Our results provide a framework for future development of cancer therapeutics targeting nutrient transport via ASCT2, as well as demonstrate the utility of combining computational modeling and cryo-EM for solute carrier ligand discovery.

An Overview of Cell-Based Assay Platforms for the Solute Carrier Family of Transporters.

The solute carrier (SLC) superfamily represents the biggest family of transporters with important roles in health and disease. Despite being attractive and druggable targets, the majority of SLCs remains understudied. One major hurdle in research on SLCs is the lack of tools, such as cell-based assays to investigate their biological role and for drug discovery. Another challenge is the disperse and anecdotal information on assay strategies that are suitable for SLCs. This review provides a comprehensive overview of state-of-the-art cellular assay technologies for SLC research and discusses relevant SLC characteristics enabling the choice of an optimal assay technology. The Innovative Medicines Initiative consortium RESOLUTE intends to accelerate research on SLCs by providing the scientific community with high-quality reagents, assay technologies and data sets, and to ultimately unlock SLCs for drug discovery.

Dual targeting of salt inducible kinases and CSF1R uncouples bone formation and bone resorption.

Bone formation and resorption are typically coupled, such that the efficacy of anabolic osteoporosis treatments may be limited by bone destruction. The multi-kinase inhibitor YKL-05-099 potently inhibits salt inducible kinases (SIKs) and may represent a promising new class of bone anabolic agents. Here, we report that YKL-05-099 increases bone formation in hypogonadal female mice without increasing bone resorption. Postnatal mice with inducible, global deletion of SIK2 and SIK3 show increased bone mass, increased bone formation, and, distinct from the effects of YKL-05-099, increased bone resorption. No cell-intrinsic role of SIKs in osteoclasts was noted. In addition to blocking SIKs, YKL-05-099 also binds and inhibits CSF1R, the receptor for the osteoclastogenic cytokine M-CSF. Modeling reveals that YKL-05-099 binds to SIK2 and CSF1R in a similar manner. Dual targeting of SIK2/3 and CSF1R induces bone formation without concomitantly increasing bone resorption and thereby may overcome limitations of most current anabolic osteoporosis therapies.

Crowdsourced mapping of unexplored target space of kinase inhibitors.

Despite decades of intensive search for compounds that modulate the activity of particular protein targets, a large proportion of the human kinome remains as yet undrugged. Effective approaches are therefore required to map the massive space of unexplored compound-kinase interactions for novel and potent activities. Here, we carry out a crowdsourced benchmarking of predictive algorithms for kinase inhibitor potencies across multiple kinase families tested on unpublished bioactivity data. We find the top-performing predictions are based on various models, including kernel learning, gradient boosting and deep learning, and their ensemble leads to a predictive accuracy exceeding that of single-dose kinase activity assays. We design experiments based on the model predictions and identify unexpected activities even for under-studied kinases, thereby accelerating experimental mapping efforts. The open-source prediction algorithms together with the bioactivities between 95 compounds and 295 kinases provide a resource for benchmarking prediction algorithms and for extending the druggable kinome.